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We study the issue of simultaneous estimation of several phase shifts induced by commuting op- 
erators on a quantum state. We derive the optimal positive operator-valued measure corresponding 
to the multiple-phase estimation. In particular, we discuss the explicit case of the optimal detection 
of double phase for a system of identical qutrits and generalise these results to optimal multiple 
phase detection for d-dimensional quantum states. 



I. INTRODUCTION 

The issue of phase estimation has important applications in quantum computation and quantum information theory. 
For example, it was shown that the existing quantum algorithms can be described in a unified way as quantum 
interference processes among different computational paths where the result of the computation is encoded in a phase 
shift [1] . The design of optimal phase measurement procedures is also crucial in various tasks of atomic physics, such 
as for example methods for precision spectroscopy [2], and quantum interferometric experiments in quantum optics. 

The problem of the optimal estimation of the value of a phase shift experienced by a quantum state has been 
extensively studied [3] , and in particular a method to derive the optimal measurement procedure in the phase covariant 
case was reported in Ref. [4]. More recently, a general formulation of the phase estimation problem was considered 
in Ref. [5], where a method to derive the optimal positive-operator valued measurement (POVM) [6] for a generally 
degenerate phase shift operator was developed. In this work we introduce the problem of simultaneous estimation of 
several phase shifts undergone by a quantum physical system. More specifically, we address the problem of estimating 
the values of M independent phase-shifts (j = 1, M), pertaining to the unitary transformation 

p^,^.^=e-S".*^-^>oe'i:".*^^^- (1) 

where Hj represent M commuting self-adjoint operators, which are in general degenerate on the Hilbert space Ti. 
of the considered quantum system and each of them has a discrete spectrum Sj {Sj can be for example Z, N, or 
Zq, g > 0) [7]. In Eq. (1) po is a generic initial pure state |V'o)(V'o I describing a quantum system with arbitrary 
dimension. We want to point out that the scenario of simultaneous estimation of several phases may be useful to 
improve the efficiency of quantum information processing tasks where several variables are encoded into phases in the 
same quantum states. 

The paper is organised as follows. In Sect. II we derive a general treatment of the multiple phase estimation 
problem, extending to the multi-phase case the approach presented in [5] for the case of single phase estimation. 
In Sect. Ill we derive the optimal POVM and the corresponding estimation fidelity for a system of N identically 
prepared "equatorial" three-dimensional systems. We want to point out that the possibility of encoding information 
in the states of three dimensional systems has been the object of several recent studies, for example in the context 
of quantum cloning [8] and quantum cryptography [9]. In Sect. IV we extend these results derived for qutrits to the 
case of quantum systems with arbitrary finite dimension d. Finally, we summarise and comment the results presented 
in Sect. V. 



II. OPTIMAL POVM FOR MULTIPLE-PHASE ESTIMATION 



In this section we derive the optimal POVM corresponding to the simultaneous estimation of several phase shifts 
experienced by a pure state IV'o), belonging to the Hilbert space TL, that undergoes the unitary transformation (1). 
Following the approach of Ref. [5] , we treat the estimation problem in the general framework of quantum estimation 
theory [6]. According to this framework, we first define a cost function C'({(^j}, which depends on the set of 

the M estimated values {i^j}, that are the results of the estimation procedure, and on the set of the M true values 
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{cf)j}. This function weights the errors for the estimates given the true values {<f>j}. The estimation problem 
then consists in minimizing the average cost C of the procedure defined as 

/>27r />27r /'27r /■27r p2ir p27r 

C= dcl>J d<t^2.. dcPMPoiWj}) #1/ #2../ d^MC{{^i},{'l>j})P{Wj}\{'Pj}) ^ (2) 

Jo Jo Jo JO Jo Jo 

where Po({'/>j}) is the a priori probability density for the true values and p{{4'j}\Wj}) is the conditional prob- 
ability of estimating the set of values {cfij} given the true values The average cost is minimized by optimizing 
the POVM dud^j}) which appears in the definition of the conditional probability as follows 

p{{^j}\{cl}j})d^id^2.:d^M = Tr[rf/i({^j})e"^S" 1 ^'"^ p^e'E-U <^^^^] . (3) 

In this work we consider the general scenario where all the values are a priori uniformly distributed, i.e. the 
probability density is simply given 

bypo({0i}) = (l/27r)^. Moreover, we consider the case where the errors in the 
estimates are weighted independently of the values of the phases, but they depend only on the values of (j)j — 
so that the cost function becomes an even function of M variables, i.e. C({^j}, {(/>j}) = C{{<j)j — 4>j}). From these 
requirements it follows that also the conditional probability corresponding to the optimal estimation procedure will 
depend only on the variables 4'j — 4>ji and therefore the optimal POVM will be phase-covariant, i.e. of the form 

d,m) = e-S.^.^^^^xe^S^^^^^-*^ f f ...^ . (4) 

In the above equation x is a positive operator satisfying the completeness constraints needed for the normalization of 
the POVM / dii{{(j)j}) = 1, where 1 denotes the identity operator. Actually, using Eq. (3) and the invariance of the 
trace under cyclic permutations it can be shown that p{{4>j}\{4>j}) = p{{<t>3 — 0i}) if and only if is covariant. 

Therefore the optimization of the phase estimation procedure can be performed by finding the positive operator x 
that minimises the average cost for a given cost function C{{(j)j}) and a generic initial state po- 

We will now show explicitly how to derive the optimal POVM for a broad class of cost functions and initial states 
Pq. First of all we will choose the representation where all the operators Hj are diagonal. We have assumed that 
the operators Hj commute, so we can identify a common basis of eigenvectors. The operators Hj are generally 
degenerate, and we will denote by |{nj})y a choice of (normalized) eigenvectors corresponding to eigenvalue Uj for 
the operator Hj, by ^{uj} the projector onto the corresponding degenerate eigenspace and by a degeneracy index, 
whose maximum value corresponds to the dimension of the degenerate eigenspace. 

We will now generalise the projection method developed in [5] and define Tiy as the Hilbert space spanned by 
the (normalized) vectors |{%}) oc n{„^.}|^o) 7^ with the choice of the arbitrary phases such that {{nj}\i}o) > 0. 
We can then write the Hilbert space of the system as Ti. = Ti-w (Xi 'H±, where the component Ti.^ is spanned by 
states that are orthogonal to \tpo)- Hence the POVM can be chosen of the block diagonal form on (S) ?i±, i.e. 
dfj.{{(pj}) = {{<pj}) ® d/i_L({^j}). In this way the component dfj.±{{(pj}) of the POVM acting on Wj_ can be chosen 
arbitrarily because it does not contribute to the average cost. Therefore, the optimisation of the estimation procedure 
can be performed by optimising only the component diJ,\^{{(pj}) of the POVM. 

In order to optimise the POVM we can assume that Il^nj}\'4'o) ^or all values of {uj}, since the resulting POVM 
will be optimal also for states having zero projection for some of these values. Due to the covariance property (4) and 
to the argument followed above, we can also write X = X|| ® X±- Thus, the problem reduces to finding the positive 
operator x\\ that minimizes the cost C in Eq. (2). In order to accomplish this task, we first rewrite Eq. (2) more 
conveniently as 

We will now express the operator X|| on the |{nj}) basis as follows 

X|| = l{"i})(Wi}IX{n,}{m,} ■ (6) 

{nj},{mj} 

The positivity condition for the operator x implies the inequalities 

|X{nj}{mj}| ^ ^^X{nj}{nj}X{mj}{mj} 1 ; (7) 



2 



where the last equahty comes from the POVM completeness relation J d^\\{(f)) = 

The cost functions we will consider are 27r-periodic functions in the variables and therefore they can be written 

as 

Ii,12,---Im=—oo 

with the condition c^i.y = c^-i^} due to the fact that the cost is a real and even function. By performing the integrals 
in Eq. (5) and exploiting the relation J^^ d^e*^""™)"^ = S„^rn/'2'^ we arrive at the following form of the average cost 

C = -Co- E Hh} E {M{nj}){{mj}\ipo)x{n,}{m,} , (9) 

{lj}^0 {mj-n.i}={lj} 

where the expression {Ij} ^ under the first summation symbol means that the sum contains all the values of the 
indexes Ij apart from the case where they are all zero (this contribution corresponds to the term co in Eq- (9)) and 
the expression {mj —rij) = {/j} under the second summation means that the equality — rij = Ij must hold for all 
values of the index j. 

Let us now consider the following inequality 

sign(c{i.}) E (V'o|{%-})({mj}|V'o)X{„,}{m,} < E KV'o|{nj})||({mj}|V'o)| , (10) 

{mj-nj}={lj} {'mj-nj}={lj} 

where we remind that we chose ({% }|V'o) > 0. The above relation becomes an equality if the conditions 

X{n,}{m,} = sign(c{„._„.}) (11) 
can be fulfilled. In this case the minimum cost takes the simple form 

C'=-co- E ^ft} E KV'olKDIKim.ll^o)! . (12) 

Notice however that the positivity of x\\ is not generally guaranteed for any set of values of sign(c{j^}). 
Let us now define a general class of cost functions, that extends the one considered by Holevo [4], with 

Cft}>0,V{Z,}^0. (13) 

For this class the conditions (11) are trivially satisfied for all the sets of values {tij} and the optimal POVM takes 
the explicit form 

d^^\m}) = ^-^|e({0,}))(e({^,})| , (14) 
where the vectors |e({^j})) are defined as 

K{</'i})) = E^'^^"''^'l{"^})- (15) 
{"3} 

In the following sections we will illustrate more explicitly the results presented here by considering specific examples. 

III. DOUBLE PHASE ESTIMATION FOR QUTRITS 

As a simple application of the concepts presented above, let us consider the optimal double phase estimation for N 
identical three-dimensional quantum systems (qutrits) all in the state 

|V(</.,^)) = -^(|0)+e^'^|l) + e'^|2)), (16) 

where {|0), |1), |2)} represents a basis for the qutrit. In the language of the previous section we identify (j)x = (f>, (j)2 = 0, 
Hi = |1)(1 1, H2 = |2)(2 I and |^o) = (|0) + |1) + |2))/y3 for each qutrit. For the composite system of N qutrits we 
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have Hi = X^feLi U H2 = X^feLi |2)(2 where |fe denotes the projection operator onto the state \j) of 
the fc-th qutrit. The operators H\ and H2 commute, and they are diagonahsed in the basis \N — n\—n2, ni, 712)1/ of the 
states where N — m—n^. qutrits are in the state |0), ni in the state |1) and n2 in the state |2). The symbol v represents 
the degeneracy index of the corresponding subspacc, and in particular it ranges from 1 to N\/{N — ni — n2)!ni!n2!. 
Since the state of the N qutrits is symmetric under any permutation performed on the qutrits, the states defined 
in the previous section by the projection method correspond in this case to the symmetric normalised states of the 
N qutrits, which we will simply denote as |7V — ni — n2,ni,n2)s (such a state is an equally weighted superposition 
of N\/{N — Til — n2)!ni!n2! components corresponding to all the possible permutations of states with N — n\ — n2 
qutrits in the state |0), n\ in the state |1) and n2 in the state |2)). 

In this case the optimal POVM (14) for the cost functions of the generalised Holevo form (13) takes the form 

dM(0,0)^gg|e(</.,0))(e(0,0)| , (17) 

where |e(</., 0)) = E,'!,=o E!1=o^ e*("i^+"^'')|iV - m - na, ni, n2).. 

Let us now compute the fidelity of the optimal double phase estimation procedure. As a cost function we can choose 
for example 1 — -F, where F is the fidelity of the estimated state |^/'(0, 9)) with respect to the true state |V'(0) 9)). This 
cost belongs to the class (13), and therefore the corresponding optimal POVM is the one written above. This choice 
of the cost function is particularly interesting because the fidelity is the figure of merit usually adopted to describe 
other processes in quantum information theory, such as for instance cloning transformations, and therefore it allows 
a direct comparison of the efficiency of optimal phase estimation with other procedures. 

By the covariance of the procedure we can write the fidelity as 

F{(j), ^) = I (V'o ^)) 1^ = ^ [3 + 2 cos (A + 2 cos ^ + 2 cos(<^ - V)] , (18) 

where 

. N N-j 

IVo) = 7=^ E E VM{N,j, k)\N - j - k,j, k), . (19) 

In the above equation we have defined M{N,j, k) = (]^_j^k)]jm • 
The average fidelity F of the procedure is then given by 

F = j F{cP,mpodl^{ct>,n 

^iC^ C £ [3 + 2 COS + 2 cos ^ + 2 cos(.^ - ^)] | (Vo e))\^ . (20) 



By performing the integration in Eq. (20) we have 

JV N-j N-p 

3 ' 3^+2 

X [^k,q{Sj,p+l + Sj+l^p) + dj^p{dk,q+l + 6k+l,q) + dj+l^pdk,q+l + 6j^p+ldk+l,q] 



F =1 + 7^11 E E VM(Ar,j,fc)M(7V,p,g) 

j,p=0 k=a q=0 



N-iN-j-1 j— : - 

= ^ + ^E E ^Wi'^)V^^fTr^- (21) 

j=o k=o V 

We want to point out that the fidelity (21) corresponding to the optimal double phase estimation for qutrits is 
always smaller than the one for equatorial qubits, given in Ref. [10], where a single phase is estimated. 

We want also to stress that, as in the case of qubits, there is a relation between optimal double phase estimation 
and optimal cloning for states of the form (16). Actually, the fidelity of the optimal double phase estimation (21) 
for a single qutrit {N — 1) coincides with the cloning fidelity for the optimal 1 M cloning transformations, that 
take a single input equatorial qutrit and produce M output copies [11], in the limit of an infinite number of output 
copies, i.e. M 00. Moreover, this result is consistent also with the relation between optimal state estimation [12] 
and optimal cloning [13] for input qutrits whose state is completely unknown (not restricted to the form (16)). 
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Wc want to point out that other figure of merits could be considered in order to evahiate the efficiency of the phase 
estimation procedure. For example, a mean periodic "variance" V{4>, 9) = 2(sin2 (j)/2 + sin^ 9/2) could be considered 
as a cost function. In this case the cost function still belongs to the class (13), and therefore, by explicitly calculating 
the average variance in a similar way as for the average fidelity, we arrive at the form 

„ 2 N-lN-j~l I _ . _ 

/ y(</.,V)TrMM<^,^)] =2 - E M(iV,j,fc)J ^ . (22) 



IV. MULTIPLE PHASE ESTIMATION FOR SYSTEMS WITH ARBITRARY DIMENSION 

In this section we generalise the results derived above for qutrits to the case of multiple; phase estimation for systems 
with arbitrary finite dimension d (qudits). We will consider the optimal multiple phase estimation for A'' identical 
d-dimensional quantum systems all in the state 

mcPj})) = -^(|0) + e'-^^ll) + e"l>^\2) + ... + e"t'^-^\d - 1)) , (23) 

where {|0), |1), |2)...|d — 1)} represents a basis for each system. 

In the language of section II, we are considering the estimation problem for d ~ 1 phases corresponding to the 
operators Hj = |, j = 1, — 1, and j-i/'o) = (|0) + |1) + |2) + ... + \d — l))/\/d for each system. For the 
composite system of N qudits we have Hj = where as in the previous section \k denotes the 

projection operator onto the state |j) of the fc-th qudit. The operators Hj commute, and they are diagonalised in 
the basis |no,ni,n2, ...nd-i)v of the states where ng qudits are in the state |0), ni in the state |1) and so on, with 
Ej=o "'J ~ ^- ^^^^ qutrits, v represents the degeneracy index of the corresponding subspace, and in 

particular it ranges from 1 to the multinomial N\/{N — m — n2 — ...nd_i)!ni!n2!...nd_i!. Analogously to the case of 
qutrits, the POVM is optimised by choosing the symmetric normalised states of the N qudits, which we will simply 
denote as \no,ni,n2, ...124-1) g. 

The optimal POVM for the cost functions of the generalised Holevo form (13) is given by (14), with M = d—1 and 

\e{{<t>j})) = E e'^^'=^'"'*ino,ni,n2,...nd_i), . (24) 

{rij} 

In the above equation the sum over {rij} means that the variables rij take all the possible non negative values 
compatible with the constraint J2jZo ~ 

Let us now compute the fidelity of the optimal multiple phase estimation procedure derived above. As in the case of 
qutrits we choose a cost function of the form 1 — F, where F is the fidelity of the estimated state \e{{4>j})) with respect 
to the true state |e({^j})). This cost belongs to the class (13), and therefore the corresponding optimal POVM is the 
one mentioned above. By the covariance of the procedure we can write the fidelity as 

d-1 

FiWj}) = l(V'o|V({<^,}))|' = ^[d + 2^cos0,- +2 5]cos(<^,- -<^fe)] (25) 

j=l j>k 

where 



1 / 

IV'o) = V \ — i — j — j r|no,ni,n2,...nd-i)s • (26) 



VdF ri V no\ni\n2l..nd-i 
The average fidelity F of the procedure is now given by 



By performing the integrations in Eq. (27) we have 
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n n n<j-l=0 



N -ni~n2- ...nd-i , . 

ni + l ' ^ ' 

Notice that the fidehty decreases as a function of the dimension d. For example, in the case of multiple phase 
estimation on the state of a single qudit it takes the simple form 

We also want to point out that the above fidelity is larger than the fidelity of estimation of a single qudit in a 
completely unknown pure state (not restricted to be of the form (23)). Actually, the fidelity of such a universal 
procedure, which we will call Funiv,i, is given by [12] 

Funiv,l = ^_^2 ' ^"^^^ 



V. CONCLUSIONS 



In this paper we have addressed the problem of simultaneous estimation of several phase shifts induced by a unitary 
transformation acting on a quantum system. Wc have derived in a general way the optimal estimation procedure 
for an arbitrary number of phase shifts and for a wide class of cost functions. We have then specialised the results 
obtained to the case of "equatorial" qutrits and then generalised them to the case of quantum system with arbitrary 
finite dimension. An interesting result that was emphasised in this work is the connection between optimal double 
phase estimation and optimal double phase covariant cloning for qutrits. We expect also that such a connection is 
valid in arbitrary finite dimension. 

Before closing the paper, we want to point out that the scenario of simultaneous estimation of several phases 
considered in this work may be exploited to design schemes where several variables are encoded into phases in the 
same quantum states and in this way the efficiency of quantum information processing tasks may be improved. 
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